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The validity of MOND and TeVeS models of modified gravity has been recently tested by using 
lensing techniques, with the conclusion that a non-trivial component in the form of dark matter is 
needed in order to match the observations. In this work those analyses are extended by comparing 
lensing to stellar masses for a sample of nine strong gravitational lenses that probe galactic scales. 
The sample is extracted from a recent work that presents the mass profile out to a few effective radii, 
therefore reaching into regions that are dominated by dark matter in the standard (general relativity) 
scenario. A range of interpolating functions are explored to test the validity of MOND/TeVeS in 
these systems. Out of the nine systems, there are five robust candidates with a significant excess 
(higher that 50%) of lensing mass with respect to stellar mass, irrespective of the stellar initial mass 
function. One of these lenses (Q0957) is located at the centre of a galactic cluster. This system 
might be accommodated in MOND/TeVeS via the addition of a hot component, like a 2 eV neutrino, 
that contribute over cluster scales. However, the other four robust candidates (LBQS1009, HE1104, 
B1600, HE2149) are located in field/group regions, so that a cold component (CDM) would be 
required even within the MOND/TeVeS framework. Our results therefore do not support recent 
claims that these alternative scenarios to CDM can survive astrophysical data. 

PACS numbers: 95.35. +d 04.50.Kd 98.62.Sb 



I. INTRODUCTION 

Over the past few years, it has been possible to test in detail the standard ACDM model of cosmology to new levels 
of accuracy, starting the so-called era of "precision cosmology" . Hitherto the model has been very successful at fitting 
observations (see e.g. [1]). This paradigm has as its foundation two main aspects: the application of classical general 
relativity in a homogeneous Friedmann-Lemaitre-Robertson- Walker metric with positive cosmological constant A, 
along with the presence of a cold dark matter component (CDM). However, the nature of the cosmological constant 
and the type of dark matter species required are presently unknown and have been the focus of vigorous investigation 
for decades. 

In the absence of any direct detection of CDM or a fundamental theory of A, proposals have been put forward 
which go beyond the framework of ACDM and offer possible alternative approaches to fitting observations. As an 
alternative to a CDM component - first posited in order to explain the flat rotation curves of galaxies [2] - Milgrom [3] 
proposed MOdified Newtonian Dynamics (MOND) which used a simple acceleration scale modification to gravity to 
account for the unexpectedly high velocities without invoking any dark matter. MOND was constructed such that 
below a certain acceleration scale, ao m 1.2 x 10~ 10 ms~ 2 , determined by the data, the usual Newtonian gravitational 
relation for the acceleration and potential is altered to f(\a\/ao)a = — V$n, where $n is the Newtonian potential. 
The function f(x) is constrained to be positive, smooth and monotonic and it is used to interpolate between two 
gravitational regimes; standard gravity above ao (where f(x) « 1) and a stronger relation below ao (where f(x) ~ x). 

In the arena for which MOND was constructed, namely galactic rotation curves, MOND has been successfully 
applied and shown to fit the data well, though this is less clear in the case of galactic clusters. MOND was also 
provided with a relativistic partner known as TeVeS [U [5], a theory that includes additional vector and scalar 
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gravitational fields to the tensor field of general relativity. This allowed the phenomenology of MOND to be tested 
against previously inaccessible data such as the cosmic microwave background [6]. 

In previous works the validity of MOND [7J and TeVeS [HI H] models of modified gravity were tested by using 
gravitational lensing techniques, with the conclusion that a non-trivial component in the form of dark matter has 
to be added to those models in order to match the observations. In particular, in Ref. [7] an analysis of a set of 
lenses from the CASTLES survey [10] was conducted using both MOND and standard gravity. The claim of MOND 
to fit observations without dark matter was tested by comparing the stellar masses calculated from the photometry 
and the required mass of the system to fit the properties of the observed lenses. A discrepancy would indicate if the 
galaxies would still require some dark matter even with MOND, hence testing the original claim of the proposal. It 
was found in Ref. [7J that there were a number of galaxies for which significant quantities of dark matter were required 
when using MOND to account for the observed lensing. Testing at the galactic scale is particularly useful in such an 
analysis, since at this scale any contribution from massive neutrinos should be negligible, unlike at the cluster scale 
where the presence of a warmer (i.e. non-CDM) component can be added to match the observations [TTj . 

In Refs. 8, 9 a fully relativistic analysis of lensing data was conducted using TeVeS, which is one attempt to cast 
MOND in a relativistic field-theory setting. In Ref. [5], a second analysis of the CASTLES lenses confirmed that 
TeVeS also required significant quantities of dark matter. However, since the modification of gravity is dependent on 
the form of f{x) and its TeVeS equivalent n(y), a parameterised form of this function was used to explore this freedom. 
The analysis in Ref. [8] showed that the most commonly used form of /i(y) did not substantially reduce the need for 
dark matter, and this dependence increased as was altered to the form which better fits the rotation curve data. 
This implies that rotation curves and lensing may require incompatible forms of the interpolating function. This 
relation was made explicit in Ref. [9] where the parameterised fi(y) from TeVeS was fitted using rotation curve data. 
The best fit form of /i(y) from rotation curves was then shown to require significant quantities of dark matter when 
applied to strong lensing data over galactic scales, and the best fit to lensing was shown to be incompatible with the 
constraints from rotation curves. Thus, it was concluded that allowing for the freedom in the interpolating function 
would not be sufficient to prevent MOND/TeVeS from requiring dark matter, as for any one form of the function a 
mass discrepancy would be observed with either rotation curves or strong lensing data. This analysis was limited to 
a one-parameter family of interpolating functions. Many-parameter cases were not considered as the introduction of 
extra free parameters was deemed to go against the original simplicity of the MOND formalism and would have, at 
present, no theoretical motivation. 

At this stage it should be mentioned that in Refs. [HI H2] the authors claim that their non-relativistic approach 
to a similar survey as in Ref. [7J led to different conclusions, namely a successful fit for MOND without dark matter. 
As shall be discussed here, although some parts of the analysis of [12 US] were valuable and have been taken on 
board in our current study, nevertheless disagreements remain with several other parts of their work, particularly the 
discussion of their results, which is found to be somewhat misleading. Moreover, and most importantly, an extended 
analysis using new data will be presented, which indicates clearly the need for the addition of dark matter in such 
analyses in order to make MOND compatible with the data, thereby contradicting their claims. 

More specifically, the purpose of the current paper is to apply the updated analysis of strong lenses, recently 
extended to larger apertures, in order to obtain more robust conclusions about the requirement of dark matter in 
MOND and TeVeS theories. New stellar mass estimates for a number of lenses from the CASTLES survey have become 
available |14j which probe regions out to several effective radii, and thus further in the deep MONDian regime and 
away from the baryon-dominated core. In the ACDM paradigm, such regions should be dominated by dark matter. 
Therefore, these data provide an opportunity to rigorously check the conclusions of Refs. [7J|H] where dark matter was 
found to be needed by MOND and TeVeS. It must be stressed here that the authors of Ref. [H] have pointed out some 
issues with the non-relativistic analysis presented in Ref. [7j. Along with the extended lensing survey discussed here, 
these issues can be addressed and their effects on the previous conclusions can be investigated in detail. As shall be 
demonstrated in this work, the use of the new expanded data and the incorporation of the suggested refinements to 
the methodology, will allow the conclusions about the confrontation between MOND/TeVeS and gravitational lensing 
observations to be made with greater certainty than had been previously possible, contradicting the conclusions and 
anticipations of Ref. [13] • 

The structure of this article is as follows: after a short review of the theory of MOND/TeVeS, outlined in section [TT| 
we focus in section III on gravitational lensing and the role of the choice of the interpolating function. A comparison, 
between lensing and stellar masses is presented in section |IV[ for a sample of strong gravitational lenses at galactic 
scales. The discussion is completed in section [V] by drawing the conclusions of this analysis. 
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II. THEORY OF MOND AND TEVES 



MOND proposes a modified relation between the Newtonian gravitational potential and the acceleration, namely 



(2.1) 



Under spherical symmetry, this relation can be rewritten in the two following forms relating the modified acceleration, 
a and the Newtonian acceleration, a^: 



a(r)f 



o, 



S(r) = f 



_ 7-1/2 
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(2.2) 



The interpolating function f(x) is more commonly used in the MOND literature. However, for the purposes of 
calculating the deflection angle of light in MOND, f(x) is more useful. These two forms of the interpolating function 
are related in the following way 



f(xf(x)) = f(x) 



(2.3) 



The above functional equation is used to convert from f(x) to f{x). 

For TeVeS, a fully relativistic analysis requires a derivation of the modified Einstein equation from the Lagrangian 
and solving the equation under the assumption of a spherically symmetric metric. For details the reader is referred 
to Refs. [H |H1 EB UH ESI ■ TeVeS [3] is a bi-metric model in which matter and radiation do not feel the Einstein metric, 
g a p, that appears in the canonical kinetic term of the (effective) action, but a modified "physical" metric, g a p, related 
to the Einstein metric by 



-20, 



where U, 



g a /3 = e "''gap - UaUfsie^ - e zcp ) , 
, - Y denote the TeVeS vector and scalar field, respectively. The TeVeS action reads 
f 



(2.4) 
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(2.5) 



where k and K are the coupling constants for the scalar, vector field, respectively; I is a free scale length related to ao 
(c.f below); a is an additional non-dynamical scalar field; J-^ v = Uu V — U u ^; A is a Lagrange multiplier implementing 
the constraint g a ^UaUp — —1, which is completely fixed by variation of the action; the function F(kGa 2 ) is chosen 
to give the correct non-relativistic MONDian limit, with G related to the Newtonian gravitational constant 120) . 
Gn, by G — [47tGn(2 — K)]/[{2 — K)k + Sir]. Covariant derivatives denoted by | are taken with respect to g M „ and 
indices are raised/lowered using the metric g^ v . 

The modified equations of motion can be calculated from the Lagrangian. For the modified Einstein equation it is 
found [1 H5] 



G a fi + 9*pA = 8ttG T al 3 + (1 - e-^)U^ {a U p) + r a p 
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For the vector field it is obtained 
87rG(f - e" 
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and similarly for the scalar field, namely 

Hy)(g ap - u a u^, a y = kG[ g a ^ + (1 + e -^)u a u^f a0 , (2.9) 

with fj,(y) defined by 

H{y) = kGa 2 , (2.10) 

where 

1 2 dF(fx) 

j/ = M 2 ^" -U»U»)cf) jfl <t> !V . (2.11) 

The function ^(y) plays the same role in TeVeS as the f(x) interpolating function does for MOND. Most choices of 
fix) can easily be converted into the TeVeS /i(y) counterparts. 

In Rcfs. [8, 9 and here, a spherically symmetric metric is assumed, which is motivated by the spherical symmetry 
of the mass profiles of the galaxy samples used in the analysis. The most general form of such a metric reads: 

g a pdx a dx fi = -e v dt 2 + e c (dr 2 + r 2 d9 2 + r 2 sin 2 9d<p 2 ) , (2.12) 

where v and £ are both functions of r. Isotropy makes the scalar field dependent only on r, namely <p — 4>{r). Matter is 
approximated as an ideal pressure-less fluid, T a p — pu a up. On the assumption that the time-like vector field has only 
one non-zero temporal component, as required by the isotropy of the Universe, the normalisation condition imposed 
by the Lagrange multiplier restricts the vector field to be 

U a = (e-" /2 , 0,0,0) . (2.13) 

In this paper we shall not consider perturbations of the TeVeS fields. Such perturbations carry important impli- 
cations over cosmological scales, as discussed for instance in Refs. [13 HH] ; where it was argued that for non trivial 
scalar TeVeS fields, depending on the cosmic scale factor, vector perturbations can play an important role in galaxy 
growth, thereby mimicking dark matter models. Our analysis probes significantly smaller (i.e. galaxy) scales, where 
we consider only local static solutions of the TeVeS fields, and thus for such configurations it is unlikely that vector 
perturbations will affect our conclusions. To be capable of doing so, such perturbations must be sufficiently strong, 
but in such a case the MONDian limit of reproducing the rotation curves of galaxies, which is the raison d'etre 
of MOND/TeVeS models, would be affected significantly. Upon the inclusion of such local perturbations, one would 
probably be forced to use different intensities for different galaxies, thereby jeopardizing the homogeneity and isotropy 
of the standard cosmology, calling into question the simplicity of such a scenario over standard ACDM. 

With the above relation for the vector fields, the physical metric can be written as 

g a pdx a dx fs = -e 9 dt 2 + e c ~(dr 2 + r 2 d(9 2 + r 2 sin 2 0dip 2 ) , (2.14) 

with the quantities v and C related to v and £ by 

v = v + 2(j> ; C = C — 20 . (2.15) 

Considering the quasi-static case, the four- velocity of the fluid, u a , is taken to be collinear with U a , and then normalise 
it with respect to the physical metric, g a p, so that u a = e^U a , leading to 

f aP = ~pe 2 *U a U . (2.16) 



Thus, the scalar field equation, Eq.(2.9), along with the isotropy constraint, leads to 
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H(y) =kGe- 2 *p, (2.17) 



where a prime denotes derivative with respect to r. Upon integration, it is obtained that, 

, _ kGm s (< r) („ +C )/ 2 
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where a scalar mass has been defined as 



m s (< r) = 4-7T 



pe 



As shown in Ref. 0], to a good approximation, the scalar mass can be considered equivalent to the "proper" mass 
contained in the same volume. Moreover, the Lagrange multiplier appearing in the vector field, Eq.(2.13l, can be 
totally determined by the vector equation, Eq.(2.8), namely 



A = 87rG(e- 



The modified Einstein equations, Eq. (2.6), for v and £ lead to the following system of differential equations: 
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In order to solve the equations numerically, the following approximation 4J is used [51 |S] : 

kGm s {< r) k 2 G 2 m 2 (< r) 
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(2.19) 



(2.20) 



(2.21) 



Thus, the system of differential equations given above can be numerically solved, and the deflection angle of light can 
be calculated for TeVeS. For MOND the deflection angle is found by assuming, as in the standard gravitational case, 
that the deflection of light will be twice the "Newtonian" deflection. In the next section, the deflection angle within 
the framework of MOND as well as TeVeS is calculated. 



III. LENSING IN MOND AND TEVES 



A. Deflection Angle 



As presented in Ref. [3T] the MOND deflection angle equation can be derived following the same method as its 
Newtonian counterpart, and is found to be 



Acp(b) 



4Gb 



-1/2 



'GM(< y/b 2 + z 2 )\ M(< y/b 2 +z 2 ) 



«o 



b 2 + z 2 } 



[b 2 + z 2 } 



213/2 



dz 



(3.1) 



A(p(b) is the deflection angle of light, b is the distance of closest approach (which is equivalent to the impact parameter 
in non-relativistic systems), z is the distance along the line of sight to the lensing galaxy with b 2 + z 2 = r 2 . Thus 
M(< \Jb 2 + z 2 ) — M(< r) and this is a cumulative mass profile. Here it has been assumed that the deflection of 
photons is twice that of non-relativisitc particles and that the photon path is nearly linear. The above relation is the 
same as that for standard gravity, except for the inclusion of the f(x) function. 



For TeVeS, the deflection angle of light is found by using the form of our metric, Eq. (2.14|, to derive the equation 
for the deflection of light in the physical metric. It reads [5] 



Aip = 2 



1 



,C(r)-S(r): 



b 2 



-1/2 



dr 



(3.2) 



where 6 is the distance of closest approach for the incoming light ray and it is related to ro, the impact parameter 
through 



'0 



(3.3) 
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The constants in the TeVeS action, Eq. (2.5), need to be fixed; they are taken to be 

k = 0.01 ; K = 0.01 ; I = Vfc£/(47rSa ) ; 
cj) c = 0.001 ; 3 = 1 + K/2 - 20 c . (3.4) 

The values of K and k are constrained [3] from solar system tests on gravity to be < 0.1, and by rotation curves to be 
> 0.001. The scale t is related to the MONDian acceleration scale, ao and b. The latter quantity is found by taking 
the limit of the function y(ii) when which then takes the form [32] y(/u) ~ bfi 2 , so for the class of fi functions 

considered here, b = 3 is set. Finally it is noted that, for <p c , the present day value of scalar field at cosmological 
scales, there are no tight constraints on its exact value, with an approximate upper bound coming from cosmological 
data. 

However, before the deflection angles of both MOND and TeVeS can be calculated, the final part which needs to 
be considered for both theories is the form of the interpolating function. 



B. The Interpolating Function and Lensing Mass Estimates 

A parameterised form of the MOND and TeVeS interpolating function will be used for the analysis presented here, 
as there is a large degeneracy in the acceptable forms that this function may take. The most basic one, first used 
with MOND, is referred to as the "simple" form. It is given by 

/w 

k V 3 

This interpolating function is not often used, as with rotation curves it does not give as good a fit as the "standard" 
MONDian interpolating function 

However, this function unfortunately becomes multi-valued when converted to the TeVeS (J>(y) interpolating function. 
Finally, the last commonly used definition is the "toy" function developed by Bekenstein in Rcf. |4_ , 

/ w° ird j 7TTS ; "M-Vl- (3 ' 7) 

The TeVeS fx(y) is not precisely the function given by Bekenstein, but approximates this function in the regime where 
our analysis takes place. However, it was noted in Refs. [531 IM] that this function gives worse fits to the rotation 
data with respect to the standard MONDian "simple" form. The authors suggested a parameterised form of the 
interpolating function, 

2x a/ 1 " 
1 + (2 - a)x+ - x) 2 + Ax l-^Vf 

This parameterised interpolating function reproduces the "simple" form when a = 1, and Bekenstein's function when 
a = 0. It was used in Ref. [9] with rotation curves, where it was found that the case a = was incompatible with 
the data (as already hinted at in Refs. [23JI21]). It was also found that for a Chabrier IMF, a = %M\ 4 2 f was the 
best fit value for rotation curves and with a Salpeter IMF a = 11.56§ 4 fg 3 gave the best fit (uncertainties quoted at 
the 95% confidence level). This range of parameters will be used as the basis of our analysis. Negative values of a 
are not considered due to their especially poor fit with rotation curves. 

The authors of Ref. [13] noted that in Ref. [7] the deflection angle relation used there made the assumption that 
f(x) = f~ 1 / 2 (x). They argued that this assumption only holds exactly for one particular form of the interpolating 
function, where f(x) = 1 for all accelerations above ao, followed by a change to f(x) = x for accelerations at or below 
ap. For other more natural forms of the interpolating function, this assumption would not give the exact result. Thus 



the proper treatment requires the use of Eq. (2.3) to convert between the two functions, and this relation will be used 



throughout the remainder of the analysis in this paper. 

Once a particular form of the interpolating function is chosen, Eqs. (3.1 ), (3.2 ) can be used to calculate the deflection 
caused by any mass profile. For the analysis presented here, a Hernquist profile [25, is applied to describe the mass 
distribution in the lenses. This profile has a cumulative mass profile of the following form 



(3.9) 
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where M is the total mass of the galaxy out to r = oo and is the core radius scale. The is related to the 
projected two dimensional effective radius by R e = 1.8153 rv The effective radius is usually defined as the projected 
aperture that contains one half of the total observed flux. For the sake of simplicity, it is assumed that there are no 
significant radial trends in the mass-to-light ratio of the underlying stellar populations, so that light can be directly 
mapped into mass. We also re-define R c as the projected radius that contains half of the stellar mass of the galaxy. 
The Hernquist profile is often used to model the baryon mass distribution of early-type galaxies, such as those which 
appear in the CASTLES survey, because its projection reproduces the characteristic R 1 / 4 surface brightness profile 
typical of such galaxies [25] . 



The lens equation allows us to analyse the deflection angle, Aip (given in Eqs. (3.1 1, (3.2)) in terms of the image 



positions 8, the source position f3, the total mass of the lens, M, and the geometry of the system, as follows 

[3 = 9- A<p(9, M, b)— . (3.10) 

Here D\ s and D s are the angular diameter distances between the lens and the source and the observer and the source, 
respectively. The actual position of the source, /?, is an unknown, as is the total lensing mass M. The image positions, 
9, are measured from the images, and the angular diameter distances are calculated from the measured redshifts of 
the lens and the source. However, the calculation of these distances is weakly dependent on the cosmological model 
assumed. Here a concordance cosmological model of (f2 m ,f2A,^k) = (0.3,0.7,0) is used and, as shown in Ref. [8], 



deviations from this model only lead to minor changes in the analysis. Using Eq. (3.10), there are two unknowns; j3 



and M. In order to solve for these unknowns, both images in the double lens systems are used. The resulting total 



mass is then used along with Eq.(3.9) to obtain a projected 2D mass within a given aperture. It is also noted that 
both stellar and lensing mass profiles can only be constrained as projected mass distributions, obtained by integrating 
the three dimensional mass densities along the line of sight. Reference [14] recently published stellar and lensing mass 
estimates out to larger radii than previously available, allowing us to extend our previous analysis (Refs. [IHS]). These 
new estimates therefore probe deeper in the regions where the baryon density is low and dark matter is dominant, 
within the ACDM framework. The study of the mass profile in strong galaxy lenses out to large apertures therefore 
imposes strong constraints on MOND/TeVeS over galaxy scales. 

In this paper, two alternative choices for the stellar initial mass function (IMF) are considered. The IMF is defined 
as the distribution of stellar masses at birth, and for our purposes it is relevant as it dominates the conversion of 
light into mass. It is usually assumed to be a universal function, although variations have been recently claimed for 
massive galaxies, where a significant excess of low-mass stars are deemed responsible for the observed strength of 
some spectral features [35], or the kinematics of the stars |27) . For this reason the analysis presents two choices of 
IMF, a classical Salpeter function [55], which consists of a single power law, therefore along the lines of the claimed 
excess of low-mass stars; and a Chabrier [29] IMF, which truncates the power law with a lognormal distribution at 
the low mass end, resulting in systematic lower values of the mass to light ratio. In the central regions of galaxies - 
where even within the standard ACDM framework no significant amounts of dark matter are expected - comparisons 
of lensing and stellar masses have been found to depend critically on the choice of IMF, with some choices of the 
mass function giving unphysically (in the absence of dark matter and assuming a pressureless fluid) higher stellar 
than lensing masses [301 131] - However, by extending our analysis to large apertures, the domain is entered where a 
single baryon contribution for any reasonable choice of IMF is not capable of explaining the lensing data, as shall be 
shown below. 



IV. ANALYSIS AND RESULTS 



Our sample of lenses is extracted from Ref. 14] , where the projected cumulative stellar and lensing mass profiles are 
given out to large radii. Estimates of the stellar mass are essential for our comparison with MOND /TeVeS predictions 
to test whether any component in addition to the baryons is needed. The lensing mass profiles - obtained within 
the standard (general relativistic) framework - are used to test whether our approximation of spherical symmetry is 
justified for a given lens. Two key criteria are applied to select our sample from Ref. |14j . Firstly, the lens has to 
be a double image system, as quad lenses are not suited to our assumption of spherical symmetry. In addition, our 
calculation of the lensing mass must be comparable to the more complex method carried out in Ref. [14], where no 
symmetry is assumed for the lenses, and no constraint is made on the functional dependence of density with radius 
(i.e. non-parametric). Hence, only double systems are selected for which the cumulative lensing mass out to a radius 
of 2R C - calculated with our approximation of spherical symmetry and a Hernquist profile, but assuming GR - is 
compatible (within 20%) with the values quoted in Ref. |14j including error bars. This criterion is important as it 
allows us to justify the assumptions set out in this study. A final sample of nine lenses were found to be suitable for 
the analysis. One lens, Q0957, has two independent double images, and so the data for these two separate pairs of 
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Lens 


M(<2R C ) 
From Ref. [14] 


M(<2R e ) 
GR 


ASS 


HS0818 




30.07 


22.6±f 8 5 2 


BRI0952 


14.76±S:f? 


9.72 




LBQS1009 




36.01 


7Q 8+ 33 - 2 


HE1104 


72.80^;g 


79.11 


-s.oti^ 


B1152 


30.431^86 


22.6 


34.6t 2 8 5 ; 8 9 


SBS1520 


41.98±l;g 


35.49 


18.31*;? 


B1600 


16.381?;^ 


18.68 


1 9 Q+ 4 - 4 


HE2149 


26.82±^ 


20.77 


2Q 1 +10 ° 


Q0957A 


151.29±g;JJ 


153.62 


i '°-4.3 


Q0957B 


i5i.29tS:g 


157.63 


_4. n+ 3 - 3 



TABLE I: Mass estimates (in units of 10 10 M Q ) for ACDM cosmology: (fim, 12a, fik) = (0.3,0.7,0). R e is the effective radius. 
M(<2R ) is the projected lensing mass within twice the effective radius, given here from two sources. In the second column 
M(<2R e ) is given from Ref. [14] where a non-parametric general approach was used assuming standard gravity. In the 
third column M(<2R C ) is stated from our analysis assuming spherical symmetry, a Hernquist profile and General Relativity. 
Comparison of these two mass estimates for the lensing mass of the galaxies shows that, for these cases, our assumptions of 
spherical symmetry and a Hernquist profile do not introduce a significant departure (> 20% of the error bars) from the more 
complex analysis of Ref. |14] , The percentage difference between the two values of M(<2R G ) is given in the column headed by 
A%. 

images are denoted by Q0957(A) and Q0957(B). It is worth noting that, in principle, only a single, robust, case is 
needed to determine whether a MOND/TeVeS framework is valid on strong lensing systems over galactic scales. 

The comparison of our GR measurements with Ref. [T3] is given in Table [i] The second column gives the lensing 
mass estimates from Ref. |14| and the third column shows our equivalent lensing mass estimates using the same 
aperture radius. The fourth column, labelled A%, shows the difference, as a percentage, between the two lensing 
mass estimates. We now turn to estimate the lensing masses of these systems within a MOND and TeVeS framework, 
for a range of parameterisations of the interpolating functions f(x) and fj,(y), respectively. 

Our results are given in Table [n] The second column shows the value of the aperture within which stellar and 
lensing masses are estimated. This analysis extends the tabulated data presented in Ref. [14] - which was given out 
to 2R e - to the outermost regions probed by the observations. For some lenses, the geometry of the images allowed 
a measurement of the lensing mass out to ~ (5 — 6)R C . The third and fourth columns give the stellar mass estimates 
assuming a Chabrier and Salpeter IMF, respectively [14] . The following columns give the lensing mass estimates for 
GR, and different parameterisations of MOND and TeVeS. The latter does not have an equivalent of the "standard" 
interpolating function, so for this parameterisation only MOND could be used. For a — 1,0,8.54 and 11.56 both 
MOND and TeVeS were used. However only for a — 1 the results for MOND and TeVeS have been stated separately. 
This is because, as the a = 1 case shows, the difference between the MOND and TeVeS mass estimates is sufficiently 
small, so that the two results can be treated as being the same. To illustrate this point, the mass estimates for the 
a = 1 case are given up to two decimal places. 

From Table [TTJ, the stellar masses are then compared to the lensing masses to look for discrepancies in the 
MOND/TeVeS case, where, by definition, no dark matter is expected over galactic scales. The results of this analysis 
are plotted in Fig. [T]and Fig. [2] for the two different IMFs used. The error bars in the figures propagate the quoted 
uncertainties from Ref. 14J. The error on the total lensing mass within a given aperture is combined in quadrature 
with the error on the stellar mass to find the combined uncertainty on our dark matter fractions, where it is assumed 
that the error on our total lensing mass is comparable to those in Ref. [14] , 

Considering the Chabrier IMF case first, Fig. [T] shows that except for BRI0952 and B1152, all other lenses require 
an additional component to explain the lensing data. This result does not depend on the choice of interpolating 
function. Notice the gradual trend towards a higher mismatch in more massive galaxies. This has been interpreted 
within the standard context of GR as an excess of dark matter in massive galaxies (see e.g. Ref [HJ[32]). Within a 
MOND/TeVeS framework, it would not be possible to explain such trend, as the physics should be solely driven by 
the baryons. However it is noted the sample in this paper is rather small, and the different apertures used for each 
galaxy makes this interpretation difficult to claim robustly. The reader is referred to Ref. [14], where the analysis is 
made, in the context of GR, for subsamples measured within the same values of R ap /R . Note that in Ref. [9], the 
a — case was found to be ruled out by fits to six galactic rotation curves. This extended the work of Ref. [22] where 
the authors showed that a = gave a poor fit to the rotation curve of one galaxy. If the a — case is assumed to be 
unfavored by rotation curve data, then, as Fig. n] shows, eight of the nine lenses show a significant excess of (dark) 
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Lens 


Rap 


M c hab 


M sal 






M LEN s(<Rap) 














GR 


MOND 






TeVeS 






(Re) 


(<Rap) 


(<R ap ) 




Standard 


a = 1 


a = 1 


a = 


a = 8.54 a 


= 11.56 


HS0818 


5.1 


12.2 


23.4 


37.6 


33.1 


28.06 


28.14 


22.1 


35.0 


35.6 


BRI0952 


2.2 


5.2 


10.7 


10.0 


8.5 


6.77 


6.79 


5.2 


9.1 


9.3 


LBQS1009 


2.7 


7.2 


13.9 


39.3 


34.5 


28.67 


28.74 


22.3 


36.6 


37.2 


HE1104 


6.5 


18.6 


37.9 


102.1 


90.4 


74.97 


75.18 


58.2 


95.1 


96.6 


B1152 


3.4 


10.3 


19.7 


26.2 


24.0 


20.51 


20.55 


16.1 


24.9 


25.1 


SBS1520 


2.7 


12.7 


23.0 


38.9 


34.6 


28.94 


29.01 


22.5 


36.4 


36.9 


B1600 


2.1 


4.0 


7.4 


19.0 


16.9 


13.98 


14.01 


10.8 


17.8 


18.0 


HE2149 


5.4 


5.3 


10.3 


26.2 


23.2 


19.72 


19.78 


15.6 


24.5 


24.8 


Q0957 (A) 


3.4 


23.4 


41.4 


176.8 


159.4 


139.30 


139.68 


112.0 


166.9 


168.9 


Q0957 (B) 


3.4 


23.4 


41.4 


182.1 


164.2 


140.33 


140.70 


110.6 


171.5 


173.7 



TABLE II: Mass estimates (in units of 1O 1O M0 units) for a concordance cosmology: (f2 m ,f2A,f2k) = (0.3,0.7,0). R ap is the 
aperture radius measured in units of the effective radius (R e ), M c hab(<Ra P ) is the Chabrier IMF stellar mass and M sa i(<R ap ) 
is the Salpeter IMF stellar mass both calculated within the aperture radius. MLENs(<Ra P ) is the lensing mass within the same 
radius. 

matter, and taking the case when a = 8.54 - which is the best fit parameter from the rotation curve analysis with a 
Chabrier IMF - five robust candidates with an excess over 50% in mass are found. 

Considering a Salpeter IMF (Fig. [2|, the need for an additional component in all galaxies is reduced, as ex- 
pected since the stellar mass estimates are greater than for a Chabrier IMF. Once again, the best case in favour of 
MOND/TeVeS with a single baryon component is a = 0, where three galaxies (BRI0952, B1152, HS0818) would not 
require any dark matter. However note that for a Salpeter IMF, both BRI0952 and B1152 are compatible with no 
dark matter within GR. Nevertheless, the data overall shows a significant need for dark matter, and the correlation 
between the lensing mass and "dark matter" excess remains valid with a Salpeter IMF. Furthermore, a consistent 
set of five lenses shows an excess higher than 50%. Note that the rejection of MOND/TeVeS in these systems is 
strongest when the best fit from rotation curves (a — 11.56) is chosen. Hence, it can be concluded that irrespective 
of the choice of IMF, neither MOND nor TeVeS can adequately account for the lensing observed using the luminous 
baryonic matter alone; an important result which confirms our previous conclusions [7H9]. In this paper, the extension 
of the analysis to larger apertures enables us to obtain an excess of lensing mass with greater confidence, as the outer 
regions are probed where it is thought that dark matter dominates. 

The five robust lenses that pose a challenge are LBQS1009, HE1104, B1600, HE2149 and Q0957. The system 
Q0957 - the most massive one - is particularly interesting in this analysis. It consists of two independent image pairs, 
allowing its lensing mass to be calculated with two separate data sets. Not only are the lensing mass estimates from 
the different image pairs consistent, but this galaxy shows a particularly high dependence on dark matter. For the 
Chabrier case the dark matter dependence does not go below 77% (60% for the Salpeter case) even when the a = 
parameterisation is considered. However, it is noted that Q0957 is a cD galaxy in the centre of a cluster [33], which 
means that an additional component clustering over larger scales, such as a 2 eV neutrino |llj . could accommodate the 
results. However, the other four galaxies (LBQS1009, HE1104, B1600, HE2149) are located in average density regions 
(either in the field or in small groups) . Hence, for these systems the excess of lensing mass can only be reconciled if 
an additional, cold, component is added, i.e. CDM. 

Regarding those lenses with a low (even negative!) excess of lensing matter, it is noted that one should be cautious 
when interpreting those data, as the error bars in the excess are very large, and always compatible with no excess. 
Averaging over a general sample of galaxies without taking care with the effects of such values can be misleading. 
Such averaging appears in Ref. [13J. In that work, when considering a = (which is labeled as the "Bekenstein" 
parameterisation), the authors found that the average dark matter dependence for their sample was —9.9% using a 
Chabrier IMF. With the errors on the stellar masses, this result would be compatible with no dark matter. However, 
this average includes four galaxies with negative dark matter estimates (Q0142, BRI0952, Q1017 and SBS1520) which, 
according to their analysis, give dark matter contributions of —85%, —67%, —64% and —47%, respectively. If the 
errors on the stellar masses are accounted for, these figures are reduced to —15%, —29%, 47% and 11%, which means 
that only Q1017 and SBS1520 are at all compatible with positive dark matter estimates. Inclusion of these cases with 
a negative dark matter contribution, which is unphysical in a pressureless fluid baryon only scenario, will obviously 
skew the average. If one excludes these galaxies from the analysis, the average excess of lensing over stellar mass 
for the remaining sample is greater, w 27%, which would be a great deal more difficult to account for without the 
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FIG. 1: Plots of the aperture lensing mass of the systems versus the percentage discrepancy after comparison with the aperture 
stellar masses. The percentage difference between the stellar masses and the lensing masses are an indication of the dark matter 
required by the system. Mtot is the lensing mass estimates within the aperture radius, Mstar is the stellar mass estimate 
from [14] . Plot for Chabrier stellar masses. GR represents the lensing mass when calculated using standard general relativity, 
and MOND when using the modified gravity. 



addition of dark matter. For the Salpeter IMF case, so many of the galaxies are found to have a negative dark matter 
requirement that the average drops to —40%. This clearly shows that cases with negative dark matter requirements 
need to be dealt with separately. For different parameterisations, a similar issue is found in this work. The "simple" 
case (a = 1) shows an average mass excess of 13% for Chabrier and —10% for Salpeter. Removing the negative dark 
matter cases, this value increases to 42% and 26%, both sufficiently high to question the validity of MOND/TeVeS. 
However, it is emphasized that an averaging approach is not a valid method to test the consistency of these theories; 
a case by case investigation offers the clearest insight into the dark matter requirement of lensing systems. 
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FIG. 2: Plots of the aperture lensing mass of the systems versus the percentage discrepancy after comparison with the aperture 
stellar masses. The percentage difference between the stellar masses and the lensing masses are an indication of the dark matter 
required by the system. Mtot is the lensing mass estimates within the aperture radius, Mstar is the stellar mass estimate 
from |14j . Plot for Salpeter stellar masses. GR represents the lensing mass when calculated using standard general relativity, 
and MOND when using the modified gravity. 



Furthermore, Ref. [TJ] does not take into account that the analysis of Ref. [9] found a best fit parameter of the 
interpolating function from rotation curve data. Had the values of a = 8.54 for Chabrier and a = 11.56 for Salpeter 
IMFs been included, even with the averaging over negative mass estimates, a need for dark matter comparable to 
that in GR would have been found, a result which would, at the very least, soften their conclusions as regards the 
success of MOND. 
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V. CONCLUSIONS 



This paper extends previous analyses by the authors [?H9] using new lensing data, in order to perform tests of the 
MOND framework, as well as its fully relativistic generalisation TeVeS. The deflection angle equations were applied 
on a set of nine strong lensing systems from the CASTLES survey, to determine the lensing masses of those galaxies. 
The selection criteria by which galaxies were chosen was primarily dictated by our assumption of spherical symmetry. 
Only double lens systems and cases where our GR mass estimates were close to the standard mass estimates from 
Ref. |14j - where neither spherical symmetry nor any particular mass profile was assumed - were included. The 
analysis of Ref. [U] gives cumulative mass profiles out to a few effective radii. For some lenses, these calculations 
probe regions where the contribution from dark matter - in the standard framework of GR - dominates over the 
baryonic component. Hence, these data offer the opportunity to revisit and further develop our previous lensing work 
in both MOND and TeVeS and to draw firmer conclusions. It was found that, for a range of parameterisations of 
the MOND and TeVeS theories significant quantities of dark matter are present for at least five galaxies: LBQS1009, 
HE1104, B1600, HE2149 and Q0957, irrespective of the stellar initial mass function adopted. Furthermore, the choice 
of a = - which gives the lowest contribution from dark matter - is a case that has been robustly ruled out by 
rotation curve data [51 123). Other parameterisations showed a even greater contribution for dark matter. In addition 
a trend is found whereby more massive galaxies present a higher excess with respect to the baryonic mass, a result 
that would be difficult to reconcile with a MOND/TeVeS approach. 

Other options, such as a change of the cosmological parameters or variations in the values of the MOND/TeVeS 
constants have been shown to affect weakly the results [8]. The most obvious recourse for proponents of MOND and 
TeVeS theories would be to find a more complex parameterisation of the gravitational interpolating function which 
would explain the excess lensing mass found in this sample, giving up one of the main advantages of such theories, 
namely the simpler description of astrophysical systems. It is found that some of these lenses robustly rule out the 
predictions of MOND and TeVeS assuming a single (baryonic) mass component over galaxy scales. Indeed if the 
claim to universality by theories such as MOND and TeVeS is to be taken seriously, then even a single robust case 
of a galaxy shown to require dark matter would be, though perhaps not fatal, a strong challenge for these theories. 
Although it is noted that Q0957 lives in a high density region - so that the lensing analysis may be probing scales 
larger than a typical galaxy halo - still an additional set of four robust candidates are found (LBQS1009, HE1104, 
B1600 and HE2149) that pose a serious challenge to these theories. 

Thus, in sharp contrast with the conclusions of Ref. [13 which, as discussed above, had certain weaknesses in 
their analysis, our findings here support and extend our previous conclusions on the fate of MOND and at least the 
simplest of the TeVeS models. However, the future of modified theories of gravity does not end with MOND and 
TeVeS. In addition to multi-parametric interpolating functions, other modified theories of gravity such as Moffat's 
MOdified Gravity (MOG) [34] are in the early stages of being tested against lensing data [35] and, at present, appear 
to offer a promising new arena of investigation. Moreover, it is possible that more complicated models, derived from 
microscopic theories of quantum gravity, such as strings and branes, incorporate some form of modification of gravity 
at low energies, which may mimick some but not all aspects of Modified Gravity Theories. An example of such a case 
is the D-particle Universe, discussed in [36], in which low-energy modifications of the gravitational laws do occur as 
a consequence of non trivial interactions of neutral matter with D-partricle space-time defects. Such modifications 
exist in the presence of significant components of Cold Dark Matter, e.g. supersymmetric partners, which in any case 
exist in low energy limits of (super) string theories. The use of gravitational lensing in testing such more complicated 
models is an interesting avenue to pursue, which would provide complementary information to other astrophysical, 
cosmological and collider tests of such theories. 
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